Dynamic light scattering from colloidal fractal monolayers 
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We address experimentally the problem of how the structure of a surface monolayer determines the 
visco-elasticity of the interface. Optical microscopy and surface quasi-elastic light scattering have 
been used to characterize aggregation of CaCOs particles at the air-water interface. The structures 
formed by cluster-cluster aggregation are two dimensional fractals which grow to eventually form 
a percolating network. This process is measured through image analysis. On the same system we 
measure the dynamics of interfacial thermal fluctuations (surface ripplons) , and we discuss how the 
relaxation process is affected by the growing clusters. We show that the structures start damping 
the ripplons strongly when the two length scales are comparable. No macroscopic surface pressure 
is measured and this is in contrast to lipid, surfactant or polymer monolayers at concentrations 
corresponding to surface coverage. This observation and the difhculty in fitting the ripplon spectrum 
with traditional models suggest that a diff'erent physical mechanism might be responsible for the 
observed damping of ripplons in this system. 
PACS: 61.43.Hv and 68.10.Et 



I. INTRODUCTION 

Surface monolayers of synthetic surfactants, polymers, 
or biological molecules like lipids or proteins can dramat- 
ically affect the physical properties of fluid interfaces, in 
particular the surface tension and the elastic and bending 
moduli It is also known that surface properties can 
be modified by the presence of small solid particles H] 
At very low surface concentrations these macromolecules 
are, generally, in a gas phase. Their influence on the in- 
terface properties increases dramatically when an overlap 
concentration is achieved. It is of interest to be able to 
control the interface parameters, as these in turn deter- 
mine the rheology and stability of emulsions and foams, 
which are of technological relevance in different indus- 
tries from dairy processing to oil recovery. 

Cluster-cluster aggregation is an important class of 
growth processes. Particles aggregate into mobile clus- 
ters which further aggregate between each other to form 
larger clusters. In recent years a lot of work has been 
done to understand the kinetics of aggregation processes 
and the resulting structures. Well studied examples are 
aggregation of colloids, smoke particles, carbon black, 
etc While many experiments have been performed 
in 3D, it has been simpler to perform computer simula- 
tions in two dimensions. Furthermore, many aggregation 
phenomena of interest, such as those occurring on a sur- 
face, are intrinsically two dimensional. For these reasons 
there has been an effort to perform experiments also in 
two dimensions. Some experiments used the air-water 
interface to provide a planar space, and studied first the 
non-diffusive aggregation of wax balls (diameter of the 
order of a millimeter) [|j and later the diffusion limited 
aggregation of silica microspheres (0.3/im-diam) Q on 
the water surface. 




FIG. 1. Binarized optical micrographs showing the growth 
of fractal aggregates of CaCOa particles on a water surface, 
(a) 750s, (b) 990s, (c) 1830s,(d) 5640s after the creation of a 
clean surface. The scale bar in (a) is 200/xm; this is compa- 
rable to the surface ripplon wavelengths that are probed by 
SQELS. 

Further experiments created a two dimensional space 
through confinement between solid boundaries. For ex- 
ample polystyrene spheres between glass slides (l.l/xm 
and 4.7/im -diam) were studied under various aggregating 
conditions While there now exists a theoretical un- 
derstanding of the geometrical structures obtained from 
cluster-cluster aggregation, together with detailed exper- 
imental studies, there are few investigations of dynamical 
and rheological properties, and, as far as we know, no ex- 
periments have been performed in two dimensions |^]. 

We were motivated to study a system where the dy- 
namical interfacial parameters could be compared to the 
geometrical structural properties of the surface aggre- 
gates. A layer of colloidal particles undergoing aggrega- 
tion is an ideal choice as it can easily be probed in situ, 
and such a study would complement bulk experiments 
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undertaken with a similar motivation 



II. EXPERIMENT 

As a model system we have studied calcium carbonate 
(CaCOa) particles that form at the interface between air 
and a solution of calcium hydroxide (Ca(0H)2) as this 
reacts with dissolved carbon dioxide. We observe the ap- 
pearance of micron-sized particles, which are effectively 
confined to the interface. These then aggregate form- 
ing two dimensional structures that will be shown to be 
fractal. 
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FIG. 2. Time evolution of the concentration (the fraction 
of black pixels in a binarized image) (•) and the fractal dimen- 
sion D/ (t) (calculated with the "box counting" algorithm). 
The inset shows the same concentration vs. time data on 
a log-linear scale highlighting the existence of two separate 
regimes. 

This system has been described previously |^,?], and 
in it has been shown that the kinetics of aggregation 
is consistent with fast cluster-cluster aggregation condi- 
tions. Our experimental conditions are similar to those 
of but we have studied the system that is obtained 
from a 0.9g/l calcium hydroxide solution. At this con- 
centration, higher than those investigated in (0.1 to 
0.3g/l), the cluster-cluster aggregation proceeds until a 
percolating network is formed after about 25 minutes. 
The samples were prepared by mixing 20ml of 1.2g/l 
Ca(0H)2 solution with 10ml of water in a 10cm diam- 
eter petri dish at 23.0±0.1°C f\ Care was taken to avoid 



contamination by dust, air drafts and external vibrations 
during the experiments. The surface at the center of the 
sample was studied with two techniques. First, a time 
series of images of the reflected light was recorded with 
a 1200*1792 pixel digital camera (Kodak DC290) using 
a Zeiss Axioplan microscope, resulting in a resolution of 
0.75/im/pixel and a field of view up to 1.3mm. Second, 
under the same conditions, surface dynamic light scatter- 
ing (SQELS) measurements were performed, with an ap- 
paratus that is described in detail elsewhere 1 1 1 1 . Briefly, 
SQELS measures the time correlation of the intensity of 
light scattered by the ripplons on the liquid surface, in 
heterodyne conditions. Surface ripplons are the thermal 
fluctuations of a liquid surface, their amplitude is of the 
order of a few Angstroms, and their length scale ranges 
from the molecular to the system size. Incident light is 
provided by a 30mW HeNe laser, illuminating a region 
of 5mm diameter on the liquid surface. The spectrum 
of the scattered light has approximately a Lorentzian 
form and it can be described in terms of a frequency 
and a damping. It is possible in principle to analyze 
this spectrum in detail and determine the surface ten- 
sion, the elasticity of a surface monolayer, and other sur- 
face parameters. In practice this technique, which is de- 
scribed in detail in ^ and ||l^ , has been used extensively 
to study interfaces of very low tensions and the surface 
visco-elasticity of monolayers. With our setup the acces- 
sible range of wavevectors (q) is 150<q<500cm^^, cor- 
responding to interface roughness wavelength A — 27r/q 
between 140 and 500/xm and frequencies (on water sub- 
phase) 15 < w < lOOkHz, where to = \/jq^/p, 7 being 
the surface tension and p the density of the liquid sub- 
phase. The potential of this technique for investigating 
colloid monolayers was first shown by Earnshaw [p3 . 



III. RESULTS 

In the following we are first going to discuss the char- 
acterization of the system and then the SQELS measure- 
ments of the surface fluctuations. The digital images have 
been analyzed consistently throughout the time evolution 
using Image J and Matlab routines. The first step was to 
apply a background subtraction using the "rolling ball" 
algorithm. Then the images were thresholded at a con- 
stant intensity level. These two steps provided binary 
images with a white background and the clusters marked 
as black. 



tion, which was kept in a glass container and used on the same 

«m . , ii 1 .i- i 1 • m 111 1- 1 day. The Ca(OH)2 was Aldrich 954-%A.C.S. Reagent, 

io avoid the complexities reported m 1101, doubly distilled -j \ y 

and deionized water was used to prepare the Ca(OH)2 solu- 
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FIG. 3. Average radius of gyration. An increase in cluster 
size is seen from 510s. Note that although at 990s the (Rg) 
is only about 40^m, at the same time it can be seen from 
Fig.^(b) that the clusters are very elongated and extend to 
about 200/im on their principal axis. 

Fig. |l] shows snapshots of the aggregating clusters on 
the water surface at different times. From these images 
we directly calculated the area fraction of black pixels 
and the fractal dimension D/ with the "box counting" 
algorithm (box sizes 1,2,3,4,8,12,16,32 and 64 were 
used). The clusters, clearly visible to the human eye, 
are composed of many closely spaced pixels. For stan- 
dard cluster recognition to identify them as a connected 
structure, it was necessary to perform a box 7 dilation. 
Having labeled pixels as belonging to separate clusters, it 
was possible to calculate the average size, mass, number 
of clusters for each image. It is useful to anticipate here 
that while the single particles have a diameter (growing 
from roughly 2/im to 10/im) which is very small compared 
to the surface roughness length scales probed by SQELS, 
the aggregates grow to have a comparable size and even- 
tually a single cluster spanning more than 1.3mm is ob- 
served. 

Fig. ^ shows the evolution of the concentration (the frac- 
tion of black pixels) with time. There are two regimes 
in the concentration increase. Up to 1000s the increase 
in concentration is very fast, and new CaCOa particles 
are formed. As they become visible in the microscope 
we observe their diameters to be between 1 and 3/xm. 
After 1000s the increase is slower and we believe it to 
be due to growth of the already existing particles. This 
is consistent with the estimate, from the images, that 
the mean particle diameter is around 7/im at 870s and 
around 15/im at 5640s. 

Fig. |l|(a) shows the surface after 750s, one can see that 
aggregation is already occurring and the surface is sep- 



arated in regions of higher and lower particle concentra- 
tion. Fig. ^b) shows the surface after 870 seconds when 
the presence of clusters is very clear. 
Fig. ^ shows the average radius of gyration (Rg) on each 
image. The radius of gyration of the i-th cluster is de- 
fined as usual as Rgi = J2j i^j ^ fmY^ with N 
the number of pixels, Tj the position of the j-th pixel and 
rm the position of the center of mass. The values of Rg 
on an image are binned into intervals and the average ra- 

dius of gyration {Rg) is calculated as {Rg) = YTT^^tNjk)' 
where k identifies the bin. There is no increase in {Rg) up 
to 510s, and after this it can be seen to increase roughly 
linearly with time. This growth cannot be followed be- 
yond 1500s because of the finite field of view (1.3mni). 
From 1470s onwards we observe a structure spanning 
more than 0.9mm, see Fig. |l|(c), while from 3150s on- 
wards we observe that the image is almost entirely com- 
posed of a single cluster, such as in Fig. W 



Fig. ^ shows the pair correlation function g(r) at differ- 
ent times. Clearly, as aggregation proceeds, the position 
of the minimum in g(r) shifts to bigger r and becomes 
shallower. A minimum in g(r) is an anti-correlation, in- 
dicating that on average across the images the clusters 
have a certain size. The maximum in g(r) occurs at the 
average distance between clusters. 
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FIG. 4. Pair correlation function of black pixels. From 
bottom to top the times are 750s, 990s, 1830s, 5640s, cor- 
responding to the images in Fig. hi The curves are shifted 
for comparison and the horizontal lines are at the expected 
asymptotic values. The minima in g(r) correspond to the av- 
erage cluster radius, these values are consistent with those 
from Fig. ^. The presence of secondary maxima indicates a 
loose liquid structure. 
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FIG. 5. Ripplon frequencies uj, normalized by their val- 
ues on a clean interface tj^, as a function of time. (▼) cor- 
responds to a scattering vector q — 154.5 ± 0.2cm~^, and 
= 14.4 ± O.lkHz; (•) is g = 222.6 ± 0.5cm-\ and 
uji = 26.8 ± 0.2kHz; (a) is g = 355 ± lcm"\ and 
uji = 53.9 ± QAkHz. 

This long-ranged correlated structure appears similar 
to that observed on colloidal particles undergoing 3D 
cluster aggregation |Q and to the surface aggregation 
described by Earnshaw in a 2D colloidal system [ p^ . 

In Fig.^ we show the evolution with time of the fractal 
dimension D/. This has been calculated starting from 
510s, when the clusters are first present and no deviation 
from the expected scaling in the "box counting" method 
was observed. D^ is seen to evolve rapidly from a value 
of 1.23 to around 1.6 at 1500s and to increase slowly 
thereafter as the structures restructure. Restructuring is 
observed in the form of free branches pivoting until they 
are doubly connected to the main cluster. This was also 
observed by |^ and is the probable cause for the fractal 
dimension being higher than that predicted for simple 
cluster-cluster aggregation 

The time correlation functions G(r) obtained with 
SQELS can be fitted with the form 

G(r) = B + Acos(wT -I- 0) exp(-rr) exp(-/3VV4). (1) 

The final Gaussian term is the instrumental broadening, 
which is calibrated separately on a clean liquid surface, 
and (/) is a phase term that accounts for the deviation of 
the power spectrum from a Lorentzian form . Fitting 
the data with EqJ^ yields the ripplon frequency lo and 
damping time F^ ||l[]. These parameters describe the 
interface dynamics phenomenologically. 
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FIG. 6. Ripplon damping coefficients F normalized by 
their values on a clean interface Fi, as a function of 
time. As in Fig. ^, (t) corresponds to a scattering vec- 
tor g = 154.5 ± 0.2cm"\ and Fi = 0.68 ± Q.Q&kHz; (•) is 
g = 222.6 ± 0.5cm-\ and Vi = 0.73 ± QmkHz; (a) is 
g = 355 ± lcm-\ and Fi = 2.8 ± O.lkHz. 

No knowledge of the nature of the interface is required 
for this analysis, and the fitting procedure is very sta- 
ble. This is in fact the only meaningful analysis of the 
ripplon spectrum that is possible if one does not have 
a model that provides a dispersion equation to relate 
the ripplon spectrum to the microscopic surface moduli 
|l6| . We have measured the time evolution of lu and F 
at three different scattering angles q=155, 223 and 355 
cm~^, respectively corresponding to ripplon wavelengths 
A=405, 282 and 177/xm and frequencies u;=14.4, 26.8 and 
53.9kHz. For each wavevector we have repeated the ex- 
periment three times. In Figj^ and Fig.|^ we present the 
values of uj and F versus time, divided by their values 
on a clean water surface, averaged over the independent 
evolutions and binned in appropriate time intervals. The 
values of the frequency, lu, for q=155 and 223cm~"'^ do not 
vary significantly from those of the clean water surface. 
For q=355cm^^ there is, after 1200s, a decrease of about 
6%. The damping, F, retains its clean surface value up 
to 1000s, when it begins to increase dramatically for all 
the wavevectors considered. Up to 1500s it can be seen 
that the increase in damping is bigger for bigger wavevec- 
tor. For q=355cm~^ F has a peak at about 1500s. As 
noted above, it is not possible from this data to deter- 
mine physical parameters of the suface without further 
assumptions, but some conclusions are possible. It is only 
after 1000s that the clusters have any effect on the rip- 
plons. This means that it is not simply the presence of 
particles but their aggregation into larger structures that 
affects the wave dynamics. At 1000s the diameter of the 
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clusters calculated from the radius of gyration is about 
90/im, but the clusters already extend above the capil- 
lary wave length scale at least in the direction of their 
principal axis (see Fig.0(b)). 

With the same sample preparation, the surface pres- 
sure was measured conventionally, as a function of time 
for up to 2 hours, with a Wilhelmy filter paper plate. 
Starting with the clean liquid surface, no pressure in- 
crease was measured, even after the formation of the sur- 
face layer. This is in contrast to the comparable situation 
occurring with the absorption of surface active molecules 
on a surface, where the surface layer sustains a surface 
pressure as soon as the molecules on the surface form a 
percolating film. 

It is important to try to determine the microscopic 
physical parameters that give rise to the observed uj and 
r. This is possible in many cases, as for example when 
observing the thermal ripplons on a free liquid surface 
or in the presence of a homogeneous viscoelastic surface 
monolayer. However the monolayer in the present study 
is composed of fractal clusters separated by free liquid 
surface, and it is thus heterogeneous on length scales 
comparable to those of the ripplons. We are not aware of 
a model that takes this into account. It might be that as 
a first approximation to this condition one should simply 
consider the surface as being composed of two kinds of 
regions (type a and b) , with different microscopic param- 
eters, each region scattering light with uja, and lui,,Ti,. 
In this case if uja and LOb are close to each other, we would 
expect to observe an effective broadening of the scattered 
power spectrum in time (increase in F), as the clusters 
make the surface heterogeneous (like Fig.|(b) and (c)), 
followed by a narrowing (decrease in F) as the clusters 
grow to cover the whole surface evenly (Figj^(d)). This 
qualitative behavior of F is indeed observed, in Fig.^ for 
q=355cm^^, but this cannot be considered conclusive of 
the validity of this approximation, since a uniform layer 
can lead to similar behavior of F (see for example typ- 
ical data from polymer monolayers as a function of 
increasing concentration) . For lack of a model that takes 
the heterogeneous quality of the surface into account, we 
have tried to describe the observed behavior following the 
analysis that is appropriate if a homogeneous viscoelas- 
tic monolayer is present on the surface. It is well known 
that the dispersion relation D{uj) for waves at an air- 
liquid interface, bearing a thin viscoelastic layer, is given 
by 



D{uj) = [eq^ + iujr] {q + m)~\ x 



+ iivrj {q + m) 



puj 



[ iLUj] (m — q) 



(2) 



sion and e is the dilational modulus. Solving this equa- 
tion for D(w) = gives an expression for the complex 
wave frequency, w, as a function of the scattering vec- 
tor q. The solutions describe both dilational and trans- 
verse waves. In a light scattering experiment it is only 
the transverse waves that scatter light and their power 
spectrum Pq{i-u) is given by: 



-Im 



iijjri{m -\- q) 



eq^ 



D{uj) 



(3) 



where m = ^ q^ + Re{m) > 0, 77 is the subphase 
viscosity, p is the subphase density, 7 is the surface ten- 



It is possible to fit the correlation function data with the 
Fourier transform of Eq.^, yielding the values of 7, e and 
e' directly. This approach was first introduced by Earn- 
shaw (reviewed in jl^ ) and recently followed by ourselves 
pT[ . However in the present work we have not found the 
fit with three free parameters, nor a fit with the surface 
pressure fixed to the independently measured value, to 
give consistent results. This could be because the time 
evolution of the system does not allow for enough data to 
be acquired, or because of the inadequacy of Eq.|| in de- 
scribing the surface layer. We tried, in analogy to many 
other systems Q where the viscosity diverges approach- 
ing gelation and an elastic modulus develops only after 
gelation has occurred, to fit with only e' as free param- 
eter, constraining the surface pressure to the pure liquid 
phase value (72.2mN/m) and the real part of the dila- 
tional modulus to 0. This gave us a dilational viscosity e' 
that was initially zero and increased rapidly after about 
1100s. Until the use of Eq|| is justified, this approach 
certainly cannot be considered quantitatively correct. 



IV. CONCLUSIONS 

We were motivated to study a two dimensional system 
where the rheological properties could be related to the 
structure. We have been successful in characterizing the 
aggregation process of CaCOs particles on a surface and 
measuring the effect of growing clusters on surface rip- 
plons. The data presented in Fig.^ and Fig.^ describe the 
dynamics of the surface ripplons and show that the clus- 
ters begin to modify the wave relaxation process when 
their size is of the order of the surface wavelength. On 
the basis of this data it is also clear that in the initial 
stages of the growth, a given cluster size has a damping 
effect which is stronger for bigger wavevectors, suggest- 
ing that we are probing a scale dependent gelation, where 
percolation on the length scale of a surface wave has an 
effect on that ripplon's dynamics, while the macroscopic 
behavior of the system is still liquid-like. These are new 
observations. It remains to be explained how to relate 
the dynamical behavior to the surface layer moduli. We 
have suggested what the effect of a heterogeneous surface 
could be and we have described the difficulties in fitting 
the data with the conventional surface layer model. We 
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feel that to understand the behavior shown in Fig.g and 
Fig.^ a model is needed to specifically take into account 
the fractal nature of the aggregates and the heterogeneity 
of the surface as a whole. 
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